Response of the common mode of interferometric detectors 
to a stochastic background of massive scalar radiation 
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We compute the angular pattern and the overlap reduction functions for the geodesic and non- 
geodesic response of the common mode of two interferometers interacting with a stochastic, massive 
scalar background. We also discuss the possible overlap between common and differential modes. 
We find that the cross-correlated response of two common modes to a non-relativistic background 
may be higher than the response of two differential modes to the same background. 
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I. INTRODUCTION 

It is well known that interferometric detectors can re- 
spond to the gravitational radiation in two ways, either 
through the "differential" or the "common" mode [I, 2], 
described, respectively, by the two response tensors D l £ 
and D % 1_ , 



D l l =u l v? ±v l v\ 



(1.1) 



where u, v are unit vectors aligned along the arms of 
the interferometer. All theoretical and phenomenolog- 
ical analyses presently available in the literature, how- 
ever, are almost exclusively devoted to the response of 
the differential mode. In the case of tensor (spin-two) 
gravitational radiation there is indeed no need of consid- 
ering the common mode, which is expected to be much 
more "noisy" (from the point of view of the experimental 
efficiency), and which is "blind" to the so-colled "cross" 
polarization state of tensor gravitational radiation. 

According to present schemes of high-energy unifica- 
tion of all interactions (such as supergravity and super- 
string models), however, the graviton is expected to have 
scalar partners. As a consequence, the cosmic back- 
ground of gravitational waves of primordial (inflation- 
ary) origin could be associated to a relic background of 
scalar, possibly massive, particles [3], interacting (with 
nearly gravitational strength) with existing gravitational 
antennas. The antenna response to massive scalar ra- 
diation differs in a significant way from the response to 
tensor radiation, mainly because of the different polar- 
ization states of the scalar particles with respect to the 
gravitons. 

Previous studies on the interaction of massive scalar 
waves with interferometric detectors were all concen- 
trated on the differential mode [4-6]. The aim of this 
paper is to discuss the response of the common mode 
to a stochastic background of massive, scalar radiation. 
It will be shown, in particular, that the common mode 
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can interact with scalars as efficiently as the differential 
mode and that, in the case of non-relativistic radiation, 
the correlated signal of two detectors may be enhanced 
(in some cases) for the common mode with respect to 
the differential mode, at fixed values of all the other ob- 
servational parameters. Such an enhancement could (in 
principle) compensate the higher level of noise expected 
in the context of a data analysis based on the common 
mode response. 

The paper is organized as follows. In Sect. II we com- 
pute the so-called antenna pattern functions for the re- 
sponse of the common mode to scalar and tensor radia- 
tion, comparing them with the (known) results relative 
to the differential mode. In Sect. Ill we compute the 
overlap reduction function for the cross-correlation of the 
common modes of two interferometers, while in Sect. IV 
we compute the overlap function for the cross-correlation 
of common and differential modes. In Sect. V we discuss 
a possible enhancement of the signal-to noise ratio for 
non-relativistic spectra. The main results of this paper 
are finally summarized in Sect. VI. 



II. PATTERN FUNCTIONS FOR 
SCALAR-TENSOR SIGNALS 



It is well known that the detection of a stochastic 
background of gravitational radiation requires the cross- 
correlation of the outputs of at least two detectors [1, 2]. 
The analysis of such a correlation requires the knwledge 
of the signal, or "physical strain" h(t), induced by the 
incident radiation on the single detector. Such a signal 
is obtained through the projection, on the arms of the 
detector, of the "tidal" forces which are due to the inci- 
dent radiation, and which are described by the equation 
of geodesic deviation (see for instance [7]). 

For a test particle moving on a generic scalar-tensor 
background, the equation of geodesic deviation deviation 
has to be generalized to take into account the possible 
direct coupling of the particles to the gradients of the 
scalar background (f> (in addition to the gradients of the 



2 



metric), and is given by [8] 



D 2 rf 
Dt 2 



0. 



(2.1) 



Here rf is the infinitesimal spacelike vector connecting 
two nearby (non-geodesic) trajectories, and q is the (di- 
mensionless) scalar charge per unit of gravitational mass 
of the given test particles. Perturbing the metric back- 
ground in the weak field approximation, g^ = r/^^+Sg^, 
and taking into account the possible scalar components 
of the metric fluctuations induced by the scalar field fluc- 
tuations 5<f>, we can set, in the longitudinal gauge [9]: 



= diag (2$, hij + 2^6ij) . 



X. 



(2.2) 



Here hij is the transverse and traceless (pure tensor) part 
of metric fluctuations, <f> and are the gauge-invariant 
Bardeen potentials, and X is the gauge -invariant scalar 
field perturbation. Conventions: 500 > 0, R va — R^a M ; 
and R^ a P = d^ va 13 + r w (*T va p - ...). In the non- 
relativistic and weak field limit, the relative acceleration 
of two test particles with proper separation L k , according 
to Eq. (2.1), is then given by 



C = -L k Ml Mij = 5R looJ + qdidjX = 
= ~hij + didj® - + qdidjX, 



(2.3) 



where My is the total stress tensor which includes the 
contribution of scalar and tensor radiation (see also [10] 
for the case q = 0). 

For a stochastic background it is now convenient to 
adopt a Fourier expansion of the fluctuations. Tensor 
fluctuations are massless, and can be expanded into fre- 
quency modes hij(u,h), where h is a unit vector specify- 
ing the propagation direction; scalar fluctuations are pos- 
sibly massive, and should be expanded into momentum 
modes X(p,n), with frequency v = E(p) = (p 2 + WL 2 ) 1 / 2 , 
where m — m/27r (note that we are using "unconven- 
tional" units h — 1, i.e. h — 1/2-7T, for an easier compari- 
son with experimental variables). For a unified treatment 
we set v = E = p for massless radiation, and we expand 
tensor and scalar fields as 

h ij = o / dp d 2 n e? 3 h A {p, fi)e 2mp ^ x -^ + h.c. 
2 J-00 J fio 



X = - dp d 2 h X{p,h)e 
2 J-00 Jn 2 

+ h.c. 



2ni(pn-x — Et) 



(2.4) 



(a similar expansion holds for 'J and <E>). The tensor efj, 
A = 1,2, describes the transverse and traceless polariza- 
tion states of spin-two gravitational radiation, and the 
angular integration d 2 h extends over the full solid angle 



Sl 2 - The total tidal tensor thus becomes 

1 



Mi 



2 J-00 P In 



d 2 h (2ttEY 



2 e?jh A +5ij* 



1 - — ) ^ihj^ - q-^h.hjX 



E 2 



+ h.c. 



2-Ki{pn-x—Et) 



(2.5) 



In the absence of sources for scalar fluctuations other 
than the scalar field X we can set $ = ^ [9], and the 
scalar contribution to My arising from the "electric" 
components Riooj of the Riemann tensor can be decom- 
posed as 



m 



where 



Tij — Sij fliTlj , 



Lij — fliTlj 



(2.6) 



(2.7) 



are, respectively, the transverse and longitudinal compo- 
nents of the scalar stress tensor. We shall refer to such 
a Riemannian contribution as to the "geodesic" part of 
the scalar tidal forces, while the direct contribution aris- 
ing from the scalar field gradients , —q(p/ E) 2 LijX , will 
be referred to as the "non-geodesic" part. By setting 
= —Fij, the projection over the detector response 
tensor D± leads finally to the physical strain h = FijD± , 
at the detector position xq — const, 



h(t,x ) = 



1 



+00 



dp d n 
'n 2 



F£(n)h A (p, n) + 



+Ff °(n)*(p, n) + F^(h)X(p, n) 



h.c 



J27ti{pfl-XQ — Et) 



(2.8) 



where 



*± ~ ■K L> ± e ij> 



Fr = Dl{Tij 



T 1 

E 2 ^3)> 



(2.9) 



are the tensor and scalar part of the antenna pattern 
functions [6] . The plus and minus sign refers to the com- 
mon and differential mode, respectively, defined in Eq. 
(1.1). It should be noticed that the longitudinal scalar 
strain survives even in the massless limit, because of the 
presence of the direct coupling to the gradients of the 
scalar background, didjX. 

For a graphic illustration of the various pattern func- 
tions we may consider a convenient frame in which the 
two arms of the interferometer, represented by the unit 
vectors u and v, are coaligned with the x\ and X2 axes, 




FIG. 2: Angular pattern functions for the common mode of 
an interferometer. The left panel illustrates both the tensor 
response F^* 1 of Eq.(2.15), and the scalar response -F? s of 
Eq. (2.16). The right panel illustrates the geodesic response 
F^° for relativistic scalar radiation with m/E = 1/3. 




FIG. 1: Angular pattern functions for the differential mode 
of an interferometer. The top panel illustrates the response 
to tensor radiation, F { _ +) (left) and F ( _ x) (right), ofEq.(2.13). 
The bottom panel illustrates the response to scalar radiation 
of Eq. (2.14), describing both the geodesic and non-geodesic 
pattern (which are proportional). 



respectively, and the direction n of wave propagation is 
specified by the polar and azimuthal angles and 9: 

u= (1,0,0), v = (0,1,0), 

h= (sin cos 0, sin sin 0, cos#). (2-10) 

We also introduce, in such a frame, two unit vectors x 
and y, orthogonal to n and to each other, 

x = (sin 0, — cos 0, 0) , 

y = (cos 9 cos 0, cos 9 sin 0, — sin0). (2.11) 

In terms of these vectors, the two independent tensor 
polarization states can then parametrized as 



^ij — XiXj UiVj-, 
(X) _ ~ - , . . « 



(2-12) 



obtain the usual tensor pattern functions [1], 



F 



(+) 



F 



(x) 



-Hj 
y' (x) 



-D lJ e 
2 y 



-D lJ € 
2 « 



--(1 + cos 2 9) cos 20, 

cos 9 sin 20, (2.13) 



and the (already known) scalar pattern functions [4, 6], 



Ff° = -(ilJ sin 2 9 cos 20, F^ = qF eco . (2.14) 

The angular distribution of the various pattern functions 
is illustrated in Fig. 1. In the same way we obtain, for 
the common mode -D^, the tensor patter functions, 





- 2 ^+^ 


1 

~ 2 


F (x 


- 1 n^A^ 

- 2^+ £ y 


= 0, 



and the new scalar pattern functions, 

Fl g -- 



1 + cos 2 



X 2 

m \ . 2 
, sin i 

E 



sin 2 9, 



(2.15) 



2 -(I) sin 2 ,, 
(2.16) 



For the response of the differential mode we then easily 



which are illustrated in Fig. 2. 

As already stressed in [6] , the response of the differen- 
tial mode to non-relativistic scalar radiation tends to be 
highly suppressed because of the factor p/E <C 1. Such 
a suppression is absent, however, for the geodesic part 
of the response of the common mode which, in the cx- 
tremum limit p — * 0, m — *■ E, tends to become isotropic 
and even maximal, F^ co — > 2. The important conse- 
quences of this property will be discussed in the following 
sections. 
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III. OVERLAP BETWEEN COMMON MODES 

The rest of this paper will be devoted to discuss the 
case of massive scalar radiation. We shall consider the 
correlated outputs of two independent detectors, s, = 
hi + rii, i = 1,2, assuming that the noises are statisti- 



cally uncorrelated, 







(riihj), i =/= j, and much 



larger than the physical strains, \rii\ 3> \hi\. The cor- 
related signal, over a given observation time T, is then 
given by 



(S) = 



T/2 



T/2 



dt dt' (s 1 (t)s 2 (t'))Q(t - t') 



T/2 



T/2 



dt dt 1 {hi{t)h 2 {t'))Q(t-t'), 



(3.1) 



where Q is an appropriate filter function [1], chosen in 
such a way as to optimize the signal-to-noise ratio. The 
average brackets are referred to the stochastic proper- 
ties of the scalar backgrounds X and A stationary, 
isotropic, unpolarized (massive) background is character- 
ized, in particular, by the following stochastic conditions: 



(X(p,n))=0, 



{X(p,h)X(p',h'))=0, 



{X{ Pl n)X%p\n'))=8{p-p') 5 ^^ l -S x {\p\) 

(3.2) 

(similar relations hold for the geodesic scalar background 
ty). Here Sxdpl) is the scalar "strain density", normal- 
ized as in the case of tensor fluctuations. 

The strain density can be conveniently expressed in 
terms of the spectral energy density of the stochastic 
background, dp/d\np, starting with the (canonically nor- 
malized) energy density of the scalar fluctuations, 

Px = (t °) = ^(X 2 + (d t X) 2 + m 2 X 2 ) (3.3) 

(Mp is the Planck mass). Using the expansion (2.4), and 
the stochastic conditions (3.2), we obtain 

px=47r 2_P^ dp S x {\p\)E\p), (3.4) 

from which, in units of critical energy density p c — 
3H 2 M 2 , 



n x = 



V dp x _ n 2 \p\E^ 
Pc dp 



3H 2 



Sx(\ P \), 



(3.5) 



and 



(X(p, n)X*(p', «')> - S(p - p')6 2 (h, h>) 



SH 2 



8ir 3 \p\E 2 



(3.6) 

In the limit m — > 0,p — > E = v, we may thus recover 
standard result for a stochastic gravitational wave back- 
ground [1], modulo the additional factor 1/4 present in 



the tensor case, and due to the fact that the graviton 
spectrum £l g (p) contains, for a given S^dpl), the con- 
tribution of two polaritazion states, each one of them 
associated to two possible helicity configurations, ±2. 

The above result can now be applied to the explicit 
computation of the correlated signal (3.1). For the non- 
geodesic part of the signal, due to the X background (see 
Eq. (2.8)), we obtain in particular 



QCT2 j-T/2 ,+oo 

"+°° dp 



7-oo \P\E 2 



x (p) f d 2 n F?l F^e 2vipfl - AS 



x e ~2-rri(E+E')At _|_ e 2wi(E-E')At 



(3.7) 



where Ax = x\ — x 2 is the spatial separation vectors of 
the two detectors, At = t — t', and Q(E) is the Fourier 
transform of the filter function, 



Q(t 



/+oo 
dE Q(E)e- 27rlE ^- t 'l 
-OO 



(3.8) 



Finally, F*±(n) are the scalar pattern functions of the two 
detectors. The time integration can be easily performed 
by assuming, as in the tensor case, that the observation 
time T is much larger than the typical time intervals over 
which Q^O, and extending to ±oo the limits of the dt' 
integration. This leads to the Dirac deltas S(E + E') and 
S(E — £"). Integration over dE' and dt becomes then 
trivial, and assuming Q(E) = Q(—E) we obtain the final 
result 



(S) 



where 



NT 3H ° 



16tt 3 



dp 



jl s (p)Q(E(p))Q x (p), (3.9) 



7± S W = T7 / d 2 n F?l(n) F£(n)e 



27ripn-Ax 



(3.10) 



is the so-called overlap reduction function [11] (N is 
an arbitrary normalization factor). In the case of the 
geodesic scalar background W one exactly obtains the 
same expression, with the obvious replacements fix 

i^Ff°, and 7 ^ 7 f°. 
The scalar overlap functions 7 _ have already been 
computed in [4] and discussed in [6]. In this section we 
present the results for the scalar overlap functions 7" s , 
7 + co associated to the common modes of two interfer- 
ometers. We shall consider, for simplicity, the case of 
two co-planar detectors (which is a realistic assumption 
if the spatial separation |Ax| = d of the interferometers 
is much smaller than the Earth curvature radius). The 
results are independent of the relative orientation of the 
axes of the detectors, thanks to the rotational symme- 
try of F + with respect to the polar angle 4>. Using the 
explicit definitions of F" s and F+ eo , and using the inte- 
gral representation of the spherical Bessel functions, the 
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angular integration of Eq. (3.10) leads to 



7+ (P) = 9i»^bj 



2 3 
jo (a) + —oj2(a) 



tT(p) = f 



ji(a)+ i(l) 4j2(a) 



(3.11) 



where a = 2irpd, and jo, ji,j2 are spherical Bessel func- 
tions, 



Jo (a) = 



sin a 



Jo (a) 



cos a, 



a 

2Z + 1 

J;+i(«) = J/(") " J/-i(a), I > 1- (3-12) 

a 

It is important to notice, for the further discussion, that 
the non-geodesic overlap 7 " s goes to zero as (p/E) 4 for 
p — > 0, exactly like the overlap of two differential modes 
[6]. The geodesic overlap, instead, goes to a constant, 
7+ C °(0) = 16ir/N, like in the case of resonant, spherical 
mass detectors [121. 



IV. OVERLAP BETWEEN COMMON AND 
DIFFERENTIAL MODES 

It may be useful, for phenomenological applications, to 
consider also the overlap between the common mode of 
one interferometer and the differential mode of another 
(or even of the same) interferometer. The result of such 
an overlap is strongly dependent on the relative geometric 
arrangement of the two detectors. 

Consider for instance a scalar background geodesi- 
cally interacting with the detectors, and generating the 
"mixed" overlap function 



7f- P b) 



1 

N 



d 2 n Ffl°(n) F§Z°{h)e 2mpfl - AS . (4.1) 



In the case of two co-planar detectors such an overlap 
is vanishing, quite independently of the relative distance 
and arm orientation. The overlap is vanishing also for 
detectors lying on two parallel planes, separated by an 
arbitrary distance. If, on the contrary, the arms of the 
interferometers are not lying on parallel planes, then the 
overlap of the common and differential modes is in gen- 
eral nonzero, even in the limit in which the spatial sepa- 
ration of the two central stations tends to zero, Ax — ► 0. 

For a simple illustration of this result we may consider 
a limiting configuration in which both interferometers are 
centered at the origin of the chosen frame. The arms of 
the first interferometer, u\ and v\, are aligned along the 
Xi and X2 axes, respectively, while the arms of the second 
interferometer, ui and «2, are aligned along X2 and X3, 




> y 



FIG. 3: Exampe of geometric arrangement with non-zero 
overlap of the common and differential mode. 



respectively. The two arms v\ and U2, in particular, are 
coaligned and coincident (see Fig. 3). The computation 
of the mixed, geodesic overlap function (4.1), for such a 
configuration, leads to the result 



15 E 2 



(4.2) 



This overlap goes to zero as (p/E) 2 for p — > 0: the cor- 
responding signal induced by non-relativistic radiation is 
thus suppressed with respect to the relativistic radiation 
with p = E, but the suppression is lower than in case of 
two differential modes. 

Similar results can be obtained for the mixed (i.e., 
common-differential) overlap relative to a non-geodesic 
scalar, background, 7" s _. The result is nonzero only if 
the two interferometers are lying on non-parallel planes, 
like in Fig. 3. In the non-geodesic case we recover how- 
ever the usual, strong suppression for non-relativistic ra- 
diation, 7 " s _ ~ (p/E) 4 for p -> 0. 

It should be stressed that the above results refer to the 
ideal situation in which the arms of the interferometers 
have exactly the same angular separation (in particular, 
they have been chosen to be orthogonal, according to Eq. 
(2.10)). In practice, however, we should also consider 
the case in which the angular separations of the arms 
are different. In that case, the mixed overlap may be 
nonzero even for co-planar detectors and for vanishing 
spatial separation, Ax — > 0. 

In order to illustrate this point let us consider the (pos- 
subly realistic) experimental set up of Fig. 4, in which 
the arms of the first interferometer have angular separa- 
tion [3 in the (x\, X2) plane, while the arms of the second 
interferometer are orthogonal, in the same plane, with 
the arm 112 coincident and coaligned with u\: 

Mi = 62 = (1,0,0), £>i = (cos/3,sin/3, 0), v% = (0,1,0) 

(4.3) 

Such an experimental set-up could be simply realized, at 
least in principle, by adding a third, non-orthogonal arm 
to already existing interferometers. 

The response of the second, orthogonal interferometer 
is characterized by the scalar pattern functions already 
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FIG. 4: Exampe of coincident interferometers with different 
angular separation of the two arms, and non-zero overlap of 
the common and differential mode. 



computed in in Eqs. (2.14), (2.16). For the second in- 
terferometer we obtain, from the definition (2.9), a 0- 
dependent response, 



Ff eo (/3) 



sin 



cos 20(sin 2 (3) 



■ sin 2(f>(sm (3 cos (3) 



Ff°(fi) = 2 - (|) 2 sin 2 



1 + cos 20(cos 2 (3) 



- sin 20(sin f3 cos (3) 



(4.4) 



and similar expressions for the non-geodesic pattern func- 
tions. 

Let us now compute the mixed overlap 7^ e ° at d = 0. 
Such an overlap is identically zero, independently of (3, if 
we take the common mode of the orthogonal interferom- 
eter and the differential mode of the non-orthogonal one. 
We obtain a nonzero result in the opposite case, in which 
we overlap the common mode of the non-orthogonal in- 
terferometer with the differential mode of the orthogonal 
one: 

ir-M) = i/ n ***f?w*z° - ^ (f) W 

(4.5) 

Also in this, case, however, we recover the strong, non- 
relativistic suppression factor (p/E) A . 



V. SIGNAL-TO-NOISE RATIO FOR 
NON-RELATIVISTIC BACKGROUNDS 

To discuss the sensitivity of the common mode to a 
scalar massive background we need to complete the com- 
putation of the signal-to-noise ratio (SNR), defined by 
[1]: 



SNR = (S)/AS, (AS) 2 = (S 2 ) - (S) 5 



(5.1) 



(the brackets denote the optimized correlation over a 
time T, according to Eq. (3.1)). For uncorrelated noises, 
much larger than the physical strains, the variance AS 
turns out to coincide with that obtained in the case of 
massless tensor radiation [1], 

rri r + OO 

(AS) 2 = (S 2 ) = - dv P x (\v\)P 2 (\v\)\Q(v)\ 2 , 

^ J — OQ 

(5.2) 

where Pi (\v\) are the so-called noise-power spectra of the 
two detectors, satisfying 

(m (*H (*')) = ^ / dvP^e- 2 ^-^. (5.3) 

Using the correlated signal (3.9), and switching the 
momentum integration to the frequency domain \v\ > m 
through the standard relations \v\ = \E\ — (p 2 -l-m 2 ) 1 / 2 , 
dp = (v /p)di>, we can then express the SNR as 



(SNR) 2 = 
A 



(S) 2 ^(3H 2 \ 2 (Q,A) 2 



(AS) 2 



= T 



u\ePx(w\)p 2 (w\) 

£(v) = \/ v 2 — m 2 , 



8tt 3 ; (Q,Q) 
[9(v — rri) + 6(—i> — m) 



(5.4) 



where is the Heaviside step function, and the round 
brackets denote the following scalar product [1]: 



(A,B) = 



dv A*(v)B(v)P x (\v\)P^\v\). (5.5) 



We have omitted, for simplicity, the indices "ng" and "X" 
because the above result generally applies both to the 
geodesic and non-geodesic response of the detector. The 
maximal value of the ratio (5.4) obviously corresponds, 
like in the graviton case [1] , to the choice of a filter func- 
tion Q proportional to A. The optimized value of SNR, 
written as an integral over the positive momentum do- 
main, is finally given by (see also [6]) 



SNR = 



3NVTH? 
8^ 



3NH 2 
8tt 3 



2T 



°(A,Af/ 2 = 
dp 

o p 3 (p 2 + m ) 3 / 2 
ll(p)&(p) 



Pi ( y/p 2 + m 2 ) P 2 ( sjp 2 + m 2 ) 



1/2 



(5.6) 



It must be noticed, at this point, that the typical noise 
power spectra Pi of present detectors are characterized 
by a minimum around a relatively narrow frequency band 
Az^o, centered around vq. outside that band the noises 
significantly increase, and the ratio (5.6) becomes negli- 
gible, SNR — > 0. This implies that the detectors are only 
sensitive to spectral energy densities Q(p) falling within 
the resonant band A^o- 
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For a massive background, as already noticed in [5], 
there are in principle three possibilities. If to ^ A^o 
then the noises Pi are always outside the sensitivity band, 
because P%{v) = Pi{\Jp 2 + m 2 ) ^> Pj(Ai/ ) for all modes 
p, and there is no detectable signal. If, on the contrary, 
to <C A^o, then the sensitivity band of the spectrum 
may possibly overlap with the relativistic branch of the 
spectrum (when p ~ v ~ Av ), but the non-relativistic 
branch p < to always corresponds to a very high noise 
Pi(m) » Pj(Az/o), and to a negligible signal. A reso- 
nant response to a non-relativistic spectrum can only be 
obtained in the third case, in which to ~ A^o. In that 
case, indeed, the noises keep at their minimum for the 
whole non-relativistic branch of the spectrum, because 
Pi{v) = P4y / p 2 + m 2 ) ~ Pi(m) ~ P J (At'o) for all modes 
< p < to. It is just in that case that the signal due to 
the overlap of the common modes may show an interest- 
ing enhancement. 

In order to illustrate this possibility we may consider 
a non-relativistic spectral distribution with a peak and 
sharp cut-off at p = to. For instance [6] 



n(p) = n 



p < m, 5 > 0, 



(5.7) 



with £lo = const. The spectral index 6 is positive to avoid 
the infrared divergence of the background energy density, 
J d\np fl(p) < oo (growing spectra are typical of relic 
backgrounds produced in the context of superstring cos- 
mology and, in particular, pre-big bang inflation [3, 13]). 
The SNR integral (5.6) is always ultraviolet convergent 
because of the physical cut-off associated to the maxi- 
mal energy mode amplified by inflation, and typical of 
any relic cosmic background. In the massless case the 
SNR integral is also infrared convergent, because the 
limit p = v — > necessarily leads outside the resonant 
regime v ~ A^o, and the noises blows up to infinity. 
In the massive case, on the contrary, P%{y) — > Pi(m) = 
const when p — > 0, so that the infrared behavior of the 
SNR integral is controlled by "f(p) and fl(p). For the 
given spectrum (5.7), in particular, we are lead to the 
following integral 



I± 



f 

Jo 



dpp^Siip). 



(5.8) 



If we consider the overlap of two differential modes 
then 7_(p) <~ (p/E) 4 asp — > [6], and the above integral 
is always convergent for any growing spectrum. The same 
is true for the overlap of the non-geodesic response of 
two common modes, 7+ s (p) (see Eq. (3.11)). For the 
geodesic response of the common mode, however, we have 
7+ e °(p) — > const for p — > 0; as a consequence, the integral 
becomes 



jgeo 



Jo 



dp p 



25-3 



2(5-1) 



(5.9) 



and is always divergent for flat enough spectra, S < 1, if 
extended down to p — (a similar effect is also found for 



the overlap of the monopole modes of two spherical, res- 
onant mass detectors [12]). Like in the case of spherical 
detectors such a divergence is unphysical, however, being 
removed by the fact that the physical observation time T 
is not infinite, and is then associated to a minimum (non- 
zero) resoluble frequency interval, Av = AE > T -1 . In 
the non-relativistic limit in which v = E ~ to + p 2 /2m 
this uncertainty condition defines a minimum momentum 
scale [12] 



P > Pmin = (2TO/T) 1/2 , 



(5.10) 



to be used as the effective lower bound in the SNR inte- 
gration (modes with p < p min cannot be resolved, and are 
to be included into the constant background over which 
scalar fluctuations are propagating, without contributing 
to the signal [12]). 

For the geodesic response of two common modes, and 
for flat enough spectra, for which the momentum inte- 
gral is dominated by the lower bound p — p m in(T), the 
SNR thus acquires an anomalous dependence on the 
observation time T. Using the zero momentum limit 
7 f j SCO (0) = 16ir/N of Eq. (3.11), using the spectrum 
(5.7), and assuming that the two spectral noises P\,P2 
are nearly equal to the same constant P(m) for p rang- 
ing from to to, we can obtain the estimate (modulo 
unimportant numerical factors) 



(SNR)f° ~ ^ 



to ) TO P(m) 



(toT) 1 ^/ 2 , 5<1, 



- =~E7=\( mT ) 1/2 i 6> 15-11) 

y to J to P(mj 

For 5 — 1 there is only a logarithmic correction to the 
usual time-dependence T 1/>2 . 

The result obtained for S > 1 is also valid for the 
non-geodesic response of two common modes, and for 
the geodesic and non-geodesic response of two differen- 
tial modes [6, 11], for any 5 > 0. We can see, therefore, 
that the signal associated to the geodetic response of the 
common mode grows faster with T than in all other cases, 
provided the stochastic background is characterized by a 
sufficiently flat spectral distribution of its energy density. 
A typical observation time 10 7 sec, for instance, leads to 
a relative enhancement of lO 7 ^ -5 '/ 2 , for all spectra with 
5<1. 



VI. CONCLUSION 

In this paper we have computed the angular pattern 
functions for the response of the common mode of inter- 
fcrometric antennas to a stochastic background of mas- 
sive scalar waves. We have taken into account both the 
direct (non-geodesic) interaction of the detector with the 
gradients of the scalar field background, and the indirect 
(geodesic) interaction with the scalar part of the met- 
ric fluctuations. The obtained results have been used 
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to compute the overlap functions between the common 
modes, and also between the common and differential 
modes, of two independent detectors. 

We have found that the signal-to-noise ratio induced 
by a stochastic background of non-relativistic particles 
is in general suppressed by the factor (p/E) 4 , unless we 
consider the correlation of two common modes, geodesi- 
cally coupled to the background. In that case the over- 
lap function approaches a constant in the non-relativistic 
regime, and if the scalar spectrum is flat enough (so that 
the correlation integral becomes dominated by its lower 
limit) then the signal-to-noise ratio grows faster with the 



correlation time T than in all other cases. By extending 
the observation time it becomes possible, in principle, 
to enhance the signal with respect to the corresponding 
one produced by the differential mode, possibly compen- 
sating the higher level of noise expected to affect the 
experimental analyses of the common mode data. 

It is a pleasure to thank Carlo Ungarelli for a careful 
reading of the manuscript, and for many helpful sugges- 
tions. We are also grateful to Leonardo Angclini for his 
precious help in the preparation of the angular plots dis- 
played in this paper. 
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